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Abstract 

The total lepton asymmetry I = J2f If in our universe is only poorly constrained by theories 
and experiments. It might be orders of magnitudes larger than the observed baryon asymmetry b ~ 
O(10 -10 ), \l\/b < 0(1O 9 ). We found that the dynamics of the cosmic QCD transition changes for 
large asymmetries. Predictions for asymmetries in a single flavour If allow even larger values. We 
find that asymmetries of If < 0(1) in a single or two flavours change the relic abundance of WIMPs. 
However, large lepton and large individual lepton flavour asymmetries influence significantly the 
dynamics of the early universe. 
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1 Introduction 

> ■ 

OO | In this article we want to describe the thermodynamic evolution of the standard model particles in the 
^ ; early universe. We will focus on two main events after the electroweak transition at T ew ~ 200 GeV 
[6]: the freeze out of WIMPs between 40 GeV and 400 MeV [17] and in more detail the cosmic QCD 
transition at T QCD ~ 200 MeV pQ. 

All standard model particles can be described as a particle fluid in chemical equilibrium. This 
excellent approximation holds because all interaction rates are much larger then the Hubble rate H. 
Within the Hubble time t H = 1/H particles get created and destroyed so fast, that the chemical 
^ 1 equilibrium is always fulfilled for T > few MeV. All particle interactions conserve charge-, baryon- 
^ ■ and lepton flavour number. Since our universe is charge neutral [3], the sum of all charged particles 
has to be zero. We know further from big bang nucleosythesis or the WMAP measurements, that in 
our universe are slightly more baryons then anti baryons. The difference in their densities normalized 
with the entropy density s(T) is measured as b = O(10~ 10 ) [I]. Things are more difficult for the lepton 
flavour number. In the temperature region of our interest, 200 GeV > T > few MeV, each lepton flavour 
is conserved. Neutrino oscillations would violate this assumption, but they become sufficient at much 
lower temperatures T < T osc ~ lOMeV. The lepton flavour density is denned as If — (n/ + n Uf )/s(T) 
for / = e, fi, t. rif is the net number density of the charged lepton flavour and n Uf the one of the 
neutrinos. The latter one is unknown and it might be, that the asymmetry between neutrinos and anti 
neutrinos is orders of magnitude larger then b. From big bang nucleosythesis, the cosmic microwave 
background and the large scale structure we only know |//| < 0.2 for T < T osc [5]. This assumes the 
asymmetry to be equal in all flavours l e — i M = l T . For T > T osc this does not necessarily hold. It was 
shown that initial different flavour asymmetries equilibrate due to the neutrino oscillations [7J before big 
bang nucleosynthesis. A scenario with asymmetries for example like l e = b ~ 10~ 10 but = —l T = 0(1) 
would be possible [HI [10] and in no contradiction to experimental data [5] . 
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Leaving the lepton flavour asymmetries as free parameters, we can set up the following system of 5 
equations for the evolution of the standard model particles in the early universe: 



b = y^L (2) 

= yqifii (3) 

i 

with / = e, n or r, the baryon number hi of species i, and qi = the charge of species i. In a Friedman- 
Lemaitre-Robertson- Walker universe the net number density of a particle i and its antiparticle i with 
chemical potential /ij = — m at a temperature T is given by p] 

n,(T,^) = rE{E 2 -m 2 fl 2 (f(i)-f(i)) dE 



•"' ;;>,••-' f ' ■ - - — -L ■ - V, (4) 



27r 2 7 mi v 2 ^exp^il exp^il, 

where ^ counts the degrees of freedom of the particle and the integral over the energy E runs from the 
mass rrii to infinity. The ± depends on the statistics: — for bosons and + for fermions. For fermions 
this can be approached to 

th{T,h) T>> =^ l 9l T^ + o(^)\ (5) 

Hli 



6 T \T 

T< mi n fmT\ 3/2 . (\L\ 



si^^jexp ^j. (6) 



The chemical potentials of all particles in chemical equilibrium can be expressed by five independent, 
e.g. fj, Ue , ii v , }j, Ut , fj, e , and \i u . For details see p. For a given temperature we solved (|TJ) to (EJ) for any 
lepton flavour asymmetry. With the calculated chemical potentials we calculated the energy density of 
each particle, given by 

e i = J^r E\/E 2 - m 4 2 ( ^ ) dE. (7) 

If the cooling universe reaches T ~ (l/3)m.j the energy becomes too low to create the particle % and 
it annihilates. The effective degrees of freedom g* of the particle fluid gets reduces by the number of 
effective degrees of freedom of the annihilating particle. The definition of the g* is 

30 ^ 

9* = -Zfi2^ej- (8) 
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A more detailed discussion can be found for instance in [TJ, [2] . In the following chapter we present our 
results for temperatures around the QCD transition for the case of equal distributed flavour asymmetries 
/ = J2fh — 3//. Another interesting scenario is the influence of b = I with = —l T ^> l e ~ b on the 
degrees of freedom, which is presented in chapter 3. We conclude our results in chapter 4. 



2 The cosmic QCD transition 

The standard model of particle physics predicts a spontaneous breaking of chiral symmetry of the QCD. 
This leads to a phase transition at a certain temperature Tqcd where quarks confine to hadrons. One 
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would like to have a phase diagram to locate the different states of matter. Such a phase diagram would 
than be spanned by the temperature and the baryo chemical potential \xb- Unfortunately the calcula- 
tions for such a diagram from first principles are extraordinary difficult [121 H3] • So far a combination of 
results from perturbative calculations, heavy ion collisions and lattice simulations lead to the diagram 
in figure [TJ 
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Figure 1: QCD diagram inspired by perturbative calculations, heavy ion collisions, and lattice simula- 
tions. If the trajectory crosses the dotted line, the transition is a crossover, whereas crossing the solid 
line would indicate a first order transition. The solid line may end in a critical point at (T e ,fi e ). The 
temperature of the transition is T c ~ 200 MeV [T21 [T3]. The trajectory of the universe has to start at 
high temperatures in the quark phase and end at //&(T = 0) = rriNuc, the nucleon mass in the hadron 
phase. 



For the evolution of the early universe the cosmic QCD transition is of crucial interest for our 
understanding of the early universe. One would like to have more knowledge about how quarks confine 
to hadrons, since it sets for example the starting conditions for one of the pillars of cosmology, the 
big bang nucleosynthesis. Depending on the order of the transition there might be several today 
observable relics like quark nuggets, magnetic fields, and the modification of the primordial background 
of gravitational waves [21 [H] . If the transition was of first order, the universe might have undergone a 
second period of inflation, suggested by [15]. We want to review here our results on the effect of leptons 
on the transition [1] and show, how they change figure [TJ 

Each of the five conserved quantum numbers b, If, and q in equations (tlj to ()3]) can be associated 
with a chemical potential, [is, Hq, and fii f , which can be expressed by the chemical potentials of the 
standard model particles pQ. To get an estimate of the I and b dependence of the chemical potentials we 
assume for T 3> Tqcd up, down, charm and strange quarks and the leptons massless. Solving equations 
(UJ) to ([3]) with ([5D leads to 

fi B (T > Tqcd) = (f& - l) ^ (9) 
^(T»T QCD ) = {-\h + l) S ^l (10) 
^L f (T 3> Tqcd) = (4 6 + i )^ (U) 



If the total lepton asymmetry I is orders of magnitudes larger, it becomes the dominant factor for all 
three chemical potentials. The order of \xb is then given by the order of /. We see also, that the chemical 
potentials can get negative, depending on the sign of I. For temperatures below pion mass shell we get 
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for the baryo chemical potential the following estimates: 



Hb(T < (m w /3T)) « m - Tin 



c(T) 
2bs(T) 



:i2) 



with c(T) = 4(-^) 3//2 . Independent of / all trajectories for the baryo chemical potential are unique for 
T < (m^/3T) ~ 50 MeV. All effects are shown in our full numerical calculations in figures |3] and HI For 
these pictures we calculated (DQ) to fl3) with net particle densities (TjJ for all standard model particles 
including their physical masses. 



3 Effective degrees of freedom and WIMPs 

In this chapter we want to give a brief introduction to the influence of large individual lepton flavour 
asymmetries on the relic abundance of WIMPs. A more detailed discussion can be found in [16] . 

In the following we assume the dark matter to be a single component WIMP with no asymmetry 
between WIMPs and anti WIMPs. These particles with masses m x typically between 10 GeV and 1 
TeV decouple chemically from the standard model particles between 40 GeV > Tf > 0.4 GeV [T7] 
depending on their mass. After this freeze out the WIMPs only react kinetically with the standard 
model particles. Their relic abundance Y can be approximated by [T7] 

Y ^r 45 V /2 1 as) 

1 - 1 m x M p <<r|v|> Tfc V&T 



7T 

Mp is the Planck mass, a the total cross section and v the velocity of the annihilating particles. The 
index fo indicates the freeze out time and x/ Q = m x /Tf oc ln(l/g*). 

The relic abundance depends on the relativistic degrees of freedom g* of the standard model particle 
ensemble. Any few percent effect (jHJ) would lead to a few percent effect in the relic abundance of the 
WIMPs. To get an analytic approximation of of the If dependence of jS]) we take all particles massless 
and get 

15 r°° E 3 



9 , (T , M) _^ ?ft/ ;_^ d£ (14) 



The index F denotes fermions and B bosons. Compared to the standard case where chemical potentials 
are neglected, we have additional 

If the fermion chemical potentials Hf/T are O(10 _1 ), the effect becomes non negligible. If for example 
the asymmetry in the \i- and r neutrinos is of order unity, their chemical potentials would become that 
big. We calculated OH]) numerically for the full standard model and for a scenario where J2flf = b with 
l e = b and 1^ = l T . If these asymmetries are in the neutrinos, values of fi Ufl = fi Ur = 0(0.5) are not in 
contradiction to theories [8] E] [TO] [11] nor to experiments [HE]. Our numerical calculations of [8] for the 
whole set of standard model particles show a significant effect. 
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Figure 2: The degrees of freedom calculated numerically for the scenario l e = and = —l T . For 
— 0.1 there is a significant deviation compared to the k = scenario. For ^ = 1 we observe a more 
then 20% effect. 

4 Conclusion 

Our knowledge of the relic neutrinos and the lepton asymmetries is still very poor and the possible 
effects of large lepton (flavour) asymmetries on the early universe have been overlooked so far. We 
investigated here the effect of large I and If for T ew > T > T u - osc and found that large lepton (flavour) 
asymmetries do have an observable impact on the dynamics of the early universe. 

An asymmetry b <C \l\ < 0.01 influences significantly the dynamics of the QCD transition [lj. 
Depending on the unknown phase structure of the QCD phase diagram in the fiB — T plane and the 
position of the hypothetical critical end point, a large lepton asymmetry might result in a first order 
transition. This seems to be possible for I > 0.02. 

Large lepton asymmetries determine the order of the baryo chemical potential. In general, \xb is a 
function of both, the baryon and lepton asymmetry, fi B = p, B (T,l,b). 

However, we have shown that the phase diagram for the cosmic QCD phase transition needs at 
least six dimensions, one for each chemical potential and the temperature. The structure of these phase 
diagrams needs further investigations, for instance the inclusion of interactions, to clarify the order. 

We have also introduced in this article the idea, that the unknown neutrino flavour asymmetries 
If > 1 effect the relic abundance of the WIMP dark matter. Even if the total lepton asymmetry is 
of the same order as the baryon asymmetry I = b, but individual flavour asymmetries I4 < 0(1), we 
observed a few percent effect. 

In our discussion of this work we did not include sterile neutrinos and in particular not their os- 
cillations with active neutrinos. The inclusion would possibly weaken the constraints on the neutrino 
asymmetries in the early universe. 
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Figure 3: Phase diagram for the baryo chemical potential. On the left hand is the sign dependence of 
fiB shown. In the high temperature regime a large and positive / would result in a negative trajectory, 
while a negative I leads always to a positive /is, as indicated by fl9]). In the low temperature regime of 
the hadron phase, both cases approach the nucleon mass. On the right hand picture are the trajectories 
for different I asymmetries shown. 
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Figure 4: phase diagrams for the charge and electron type lepton flavour chemical potential. The upper 
two pictures show the trajectories of the charge chemical potential and the lower the trajectories of the 
electron flavour chemical potential. Both chemical potentials depend on the same way on the sign of 
large I. Also here the splittings of the trajectories for different I are shown. 
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